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ABSTRACT 

We present hydrodynamical simulations of major mergers of galaxies and study 
the effects of winds produced by active galactic nuclei (AGN) on interstellar gas in the 
AGN's host galaxy. We consider winds with initial velocities ~ 10, 000 km s~^ and an 
initial momentum (energy) flux of ^ t^L/c (~ 0.01 t^, L), with r^, ~ 1 — 10. The AGN 
wind sweeps up and shock heats the surrounding interstellar gas, leading to a galaxy- 
scale outflow with velocities ~ 1000 km s^^, peak mass outflow rates comparable to 
the star formation rate, and a total ejected gas mass ^ 3 x 10^ Mq. Large momen- 
tum fluxes, Tu, > 3, are required for the AGN-driven galactic outflow to suppress star 
formation and accretion in the black hole's host galaxy. Less powerful AGN winds 
(tu, < 3) still produce a modest galaxy-scale outflow, but the outflow has little global 
effect on the ambient interstellar gas. We argue that this mechanism of AGN feedback 
can plausibly produce the high velocity outflows observed in post-starburst galaxies 
and the massive molecular and atomic outflows observed in local ultra-luminous in- 
frared galaxies. Moreover, the outflows from local ultra-luminous infrared galaxies are 
inferred to have r^, ^ 10, comparable to what we find is required for AGN winds to 
regulate the growth of black holes and set the Mbh — cr relation. We conclude by 
discussing theoretical mechanisms that can lead to AGN wind mass-loading and mo- 
mentum/energy fluxes large enough to have a significant impact on galaxy formation. 



1 INTRODUCTION 



Accretion onto a central massive black hole (BH) in a 
galactic nucleus produces energy in the form of radia- 
tion, relativistic jets, and wider angle (le ss-coUimated ) non- 
relativistic {v ~ 10''kms"^) outflows (|Krolik|[l999|. The 



coupling of this energy output to gas in galaxies and in the 
intergalactic medium is believed to play an important role in 
galaxy formation, potentially regulating the growth of mas- 
sive galaxies and the thermal properties of the intracluster 
medium in galaxy groups and clusters (e.g 
1998[ [Croton et al.|2006l ). 



In this paper, we assess the influence of AGN winds 
on gas in the AGN's host galaxy using three-dimensional 
numerical simulations. Previous analytic work and one and 
two-dimensional simulations have demonstrated that AGN 
winds can in principle sweep up and drive gas out of galaxies, 
potentially explaining the Mbh — o" relation and the dearth 
of gas and ongoing star formation in massive, early-type, 
galaxies (e.g., [K ing 2005: |King et al.|2011|[Novak et al.|20TT] 
Qstriker et al.|2010|and references therein). 



Silk & Rees 



The impact of this 'feedback' on the gas in galaxies is 
particularly uncertain, both because the interstellar (ISM) 
gas is denser, and thus more difficult to affect dynamically, 
and because much of the ISM subtends a relatively mod- 
est solid angle relative to a central active galactic nucleus 
(AGN). However, analytic estimates and numerical simula- 
tions have demonstrated that if a modest fraction of the 
energy produced by accretion onto a central BH can couple 
to the surrounding gas, it can unbind the interstellar gas 
(e.g., [Silk fc Rees|1998[rDi Matteo et al.|2005[ ). The physical 
processes most likely to produce such an effect are winds 
||King||2003| [King et al!]|2011[ ), radiati on pressure ([Murray] 
et al.|2005 l, and/or Compton heating ( Sazonov et al.|2004 l 
from a central AGN. Understanding how this works in detail 
is one of the major challenges in our understanding of the 
connection between AGN physics and galaxy formation. 



Observationally, there is strong evidence that AGN in- 
deed drive powerful outflows. Broad-absorption line (BAL) 
quasars, which show blue-shifted absorption lines in the rest- 
frame ultraviolet with inferred outflow velocities ~ 10, 000 — 
40, 000 km s~^, represent over ~ 40% of quasars in infrared 



selected samples ( Dai et al. 2008 1. A similar fraction of radio- 



quiet quasars show evidence for high velocity outflows in X- 
ray absorption line spectroscopy ( [Tombesi et al.|2010| ). It is 
likely that all quasars possess such outflows but that they 
are only observed when the system is viewed modestly edge- 
on ( [Murray et al.|1995[ ). However, determining the mass-loss 
rate from spatially unresolved absorption-line observations 
is notoriously difficult given uncertainties in the radius of the 
absorbing gas. In a handful of low-ionization BAL quasars 
(in particular, FeLoBALs) this degeneracy has been broken, 
suggesting mass loss rates significantly larger than the black 
hole accretion rate ([Moe et al.[[2009| [Bautista et al.[[2010| 
Dunn et al. 2010 Faucher-Giguere et al. 20111. These ob- 
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servations trace absorbers at large distances from the BH (~ 
kpc) , in contrast to most of the high ionization U V and X- 
ray absorption seen in BAL quasars, which arises at < 1 pc. 
In addition to these weU-characterized outflows, it is possi- 
ble that quasars drive even more powerful outflows during 
phases of super-Eddington accretion (King fc Pounds|2003 



King et al. 2011 1 or during obscured phases when the AGN's 



radiation is trapped in the galactic nucleus, enhancing the 
force on the ambient gas ( DeBuhr et al.|2010 2011 1. 

To characterize the impact of AGN winds on their host 
galaxy, we carry out numerical simulations of major galaxy 
mergers with models for BH growth and feedback. Although 
it is by no means certain that all quasars are associated with 
mergers, this provides a convenient and well-posed model in 
which to study gas inflow in galactic nuclei. It also allows 
us to readily compare our results to the extensive previous 
literature on AGN feedback during mergers. 

In the next section (32} we summarize our methods, em- 
phasizing our implementation of AGN winds and our treat- 
ment of ISM cooling. We then describe our key results in 
33] We conclude in !|4] by discussing the implications of our 
results for models of AGN feedback and for galactic winds 
driven by AGN. We also compare our results to observa- 
tions of high-velocity outflows from post-starburst galaxies 
and ultra-luminous infrared galaxies and summarize the the- 
oretical processes that can produce the large mass-loading 
and momentum fluxes we find are necessary for AGN winds 
to have a substantial impact on the surrounding ISM. 



2 METHODOLOGY 

We use a modifi ed version of the TreeSPH code GADGET-2 
( Springel|2005 1 to perform simulations of equal- mass galaxy 
mergers. This version of the code includes the effective star 



formation model of Springel & Hernquist ( 2003 1 (hereafter 
SII03) and the radiation pressure AGN feedback model of 
DeBuhr et "al](|2011[) (hereafter DQM). As described below. 



we modified the code further to implement a model of feed- 
back via AGN winds, and to implement a more physical 
model of how the interstellar gas cools. 

2.1 Galaxy Models and Initial Conditions 

We carry out simulations of major mergers of two equal 
mass galaxies. We simulate only a single galaxy mass and 
focus our resources on studying the effects of feedback via 
AGN winds produced at small radii. Each galaxy model has 
a rotationally supported disk of gas and stars and a stellar 
bulge, all embedded in a halo of dark matter. Both the stellar 
and gaseous disks have an exponential radial profile with a 
scale length of 3.51 kpc. The vertical profile of the stellar 
disk is that of an isothermal sheet with a scale height of 702 
pc. The vertical structure of the gas disk is set by hydrostatic 
equilibrium. The halo and stellar bulge both have |Hernquist| 
( 1990 1 profiles with a halo virial and half mass radius of 229 



and 102 kpc, respectively (a concentration of 9.0) and with 
a bulge effective radius of 1.27 kpc. Each galaxy has a total 
dynamical mass of 1.94 x lO^'^M©, a total disk mass of 7.96 x 
1O^°M0 and a bulge of mass 2.66 x IO^^Mq; 10% of the disk 
mass is gas. We have performed an additional simulation 
with a disk gas fraction of 30 %; this run gives qualitatively 



similar results and so is not discussed in detail. The black 
holes are modeled as additional coUisionless particles with 
initial masses of IQPMq. 

For our fiducial simulations, each galaxy is formed from 
2.4 x 10^ particles: the halo has 9 x 10* dark matter particles, 
the disk has 6 x 10* stellar particles and 6 x 10* gas particles, 
and the bulge has 3 x 10* stellar particles. The gravitational 
force softening length is e = 70 pc for the disk and bulge 
particles and e = 176 pc for the halo particles. We describe 
resolution tests in § 3.3. 

For the merger simulations, the two equal mass galaxies 
are placed on a prograde orbit with an initial separation of 
142.8 kpc and an initial velocity of each galaxy of 160 km 
s~^ directed at an angle of 28° from the line connecting the 
galaxies. The corresponding orbital energy is approximately 
zero. The spin of the two galaxies are not aligned with the 
orbital angular momentum of the system; the relative angle 
between the spin directions is 41°, with one galaxy's spin 
making an angle of 10° degrees relative to the orbital angular 
momentum. 



2.2 Black Hole Accretion 

In this section we briefiy review the sub-grid accretion model 
presented in detail in DQM. The sub-grid accretion rate on 
scales smaller than our resolution (both gravitational force 
softening and SPH smoothing) is estimated with a model 
motivated by the redistribution of angular momentum in 
the gas. For gas with a surface density E, a sound speed 
Cs and a rotational angular frequency Jl, the mass accretion 
rate into the nuclear region is: 



Mvisc = 37rQE- 



(1) 



Here q is a free parameter of the model characterizing not 
only the efficiency of angular momentum transport, but also 
the amount of gas that turns into stars (on scales below 
our resolution) instead of falling into the black hole. Our 
fiducial values of a range from 0.05 — 0.15, motivated (at the 
order of magnitude level) by comparison to the simulations 
of gas infiow from ~ 0.1 — 100 pc of |Hopkins fc Quataert| 
( |2010[ ). Over this range of a, there is little dependence of 
our simulation results on a (see DQM). 

To perform the sub-grid estimate of the gas properties 
in equation (ml, we take averages of the properties of the 
SPH particles inside a region of radius Race ~ 300 pc cen- 
tered on the black hole. For reasons described in DQM we 
take this radius to be four times the gravitational force soft- 
ening, e, of the particles in the simulation. 

The BHs in our simulation are modeled as specially 
marked coUisionless particles. To minimize the unavoidable 
spurious stochastic motion of the BH particles due to grav- 
itational interactions with the stellar and gas particles, we 
assign the BH particles a (fixed) tracer mass of 2 x IO^A^q, 
which is ~ lO'^-lO* times more massive than the other parti- 
cles in the simulation. During each simulation, we also com- 
pute separately the 'real' mass of the BH that grows ac- 
cording to the accretion rate at small radii (with an initial 
mass of IO'^AIq). Details of this step were given in § 2.3.2 of 
DQM. 



© 0000 RAS, MNRAS 000, 000-000 



AGN-Driven Galactic Winds 



2.3 Radiation Pressure Feedback 

The luminosity of the black hole is taken to be a fraction, 
rj = 0.1, of the rest mass energy of the accreted material: 



i(riMt„c'^,LeddY 



(2) 



The net accretion rate Afin onto the black hole need not 
equal the viscous inflow rate Afvisc at ~ 100 pc in equa- 
tion ([l| if there are significant outflows on scales < 100 pc 
(roughly our resolution); we return to this point in § 2.4. 

In the radiation pressure feedback model introduced in 
DQM, the accretion luminosity is coupled back into the sur- 
rounding gas by adding a total force 

Prad = r—, (3) 

c 

which is shared equally by all SPH particles inside Race- The 
added force is directed radially away from the BH particle. 
Here t is a free parameter of the model, representing the op- 
tical depth to infrared (IR) radiation in the nuclear region. 
Our fiducial choice for r in this work is 20. DQM showed 
that r ~ 20 was required to reproduce the observed normal- 
ization of the A/bh — o" correlation in their simulations. 

2.4 AGN Wind Feedback 

In addition to the radiation pressure feedback described by 
equation Q, accretion onto the black hole can explicitly 
drive a wind at radii well beneath our resolution. Our treat- 
ment of such a wind is motivated in particular by observa- 
tions of BAL quasars, which have ~ 10, 000 km s~^ outflows 
launched from near the broad line region at ~ 0.1 pc (e.g., 
Murray et al. [1995 1. In our model, the AGN winds at small 
radii carry a momentum flux given by 



Pw 



L 

C 



(4) 



where r,„ is a further parameter of the model representing 
the total momentum flux in the wind. The wind is launched 
at a fixed speed v^,. Thus, the rate at which mass is added 
to the wind is given by 

L 



My, = r^ 



= TwTI Mir, 

Vw 



(5) 



where we have linked the luminosity of the black hole to the 
accretion rate as in equation (pi. In the presence of signif- 
icant AGN winds, not all of the material entering the nu- 
clear region, Af„isc, actually reaches the black hole. Instead, 
the net accretion rate into the black hole, Mi„, must be re- 



duced by the mass of the outflow (see Ostriker et al.|2010 |: 
Min = Myisc — Mw, which implies that the true black hole 
accretion rate is given by 



Mi„ = min 



M^i, 



Le 



,1 + T„r;c/«„ ' rjc^ J 

Note that equation ([6|, not equation ([l]), determines the 
AGN luminosity and thus the magnitude of the feedback in 
equations ([3| and H. 

To deposit the wind momentum, we give kicks to par- 
ticles inside Race- Particles receiving a kick are selected 
stochastically with each particle having an equal probabil- 
ity of being added to the wind. The probability of a kick is 
chosen to ensure that the time average of the mass kicked 



in a given timestep At is given by M^At. The kick is im- 
plemented by adding a velocity v^ to the current velocity of 
each selected particle. The direction of this imparted velocity 
is chosen to be more heavily weighted toward the surround- 
ing disk: if 9 is the angle between the imparted momen- 
tum and the disk normal, the probability distribution over 
X — cos 6* is given by p{x) = 3(1 — a;^)/4. For this purpose, 
the disk normal is defined to be the direction of the total 
orbital angular momentum about the black hole of all the 
gas particles inside Race- This modest equatorial bias in the 



wind direction is motivated by models of BAL quasars ( Mur- 
ray et al.|1995l. Nonetheless, this choice is not that critical: 



the efficacy of the feedback is largely determined by the out- 
flow momentum/energy flux that is directed within the solid 
angle subtended by the surrounding ISM. For our fiducial 
model, ~ 40% of the momentum flux is directed within one 
scale-height of the disk at ~ 100 pc. Models with isotropic 
kicks require slightly large values of r„ to give results similar 
to our fiducial model because less of the feedback is directed 
towards the surrounding ISM. 

Despite the similarity between equations S and Q, 
our two feedback processes affect the gas particles within 
Race in different ways. In the case of our simple radiation 
pressure model, all gas particles within Race receive an addi- 
tional acceleration set by equation ([3| and are pushed out- 
ward from the BH. As shown in DQM, this feedback process 
evacuates the central gas reservoir and lowers the BH accre- 
tion rate and the final BH mass, but the velocities of the gas 
particles are not so large as to produce significant large-scale 
galactic outfiows. In our AGN wind model, by contrast, a 
small fraction ( < 5%) of gas particles within Race receives an 
instantaneous large wind velocity of 3000 to 10, 000 km s~^. 
These particles leave Race within < 10^ years and, as we 
show below, can drive a galaxy-scale outflow via their inter- 
action with the ambient ISM. 

Our implementation of wind feedback adds two addi- 
tional input parameters to the simulation. The first, v^, de- 
scribes the launch speed of the winds while the second, t^, 
describes the total momentum fiux in the wind. Our default 
wind speed is motivated in part by observations of BAL 
quasars and theoretical models of the origin of such winds 
via line driving in the accretion disk at ~ 0.1 pc (though 
our model should not be taken as a literal implementation of 
this physics). Observed wind velocities are ~ 10,000kms~^. 
Models of line driving lead to momentum fiuxes in the wind 
of < L/c, i.e., Tw < 1 because the lines do not typically 
completely cover the continuum ( Murray et al.|1995[ ). As we 
show below, however, larger values of Tu, are required for 
AGN outfiows to have a significant effect on the gas dynam- 
ics and star formation history in galactic nuclei. 

Observationally, momentum fluxes are difficult to in- 
fer in most cases because of ambiguities in the absolute 
density/radial scale at which the absorption occurs. In a 
handful of low-ionization BALs (in particular, FeLoBALs), 
this degeneracy has been broken, implying momentum fiuxes 
~ 0.3 — 5I//c, i.e., r^ ~ 0.3 — 5 (with significant uncer- 
tainties; see |Moe et al.|[2009l [Bautista et ai1[20Tol |Dunn| 
|et al.|[2010| [Faucher-Giguere et al.||2011| ). It remains uncer- 
tain whether these values of t^ are representative of the 
entire BAL-quasar population. We take t„ ~ 5 as our fidu- 
cial value but also explore the range 1 < r^j < 10. Our 
fiducial value of Tm = 5 was chosen for three reasons that 
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will become clear later in the paper: (1) the AGN wind then 
leads to a galaxy-scale outflow that significantly influences 
the surrounding ISM dynamics, (2) the final BH mass in 
the simulations is reasonably consistent with the observed 
Mbh — o" relation for r™ ~ 5 — 10. (3) Observations of high 
speed atomic and molecular outfiows in local ULIRGs sug- 
gest r™ ~ 1 ([Feruglio et al.|2010. Chung eTal.|2011|[Rupke 



& Veilleux 2011 Sturm et al. 20111. We discuss physical 



mechanisms that can produce such powerful AGN winds in 




2.5 ISM Model and Gas Cooling 

In the subgrid ISM model of SH03, gas with densities above 
the star formation threshold approaches an effective ther- 
mal energy, Ue//, set by a balance between cooling and the 
feedback from star formation. If processes such as shocks or 
adiabatic expansion/compression cause the internal energy 
to deviate from u^ff, the differences decay on the timescale 
given by eq. (12) of SII03. This decay timescale is set by the 
subgrid model rather than the true cooling time of the gas. 
One consequence of SH03's ISM model is that suffi- 
ciently dense shock heated gas does not cool on its cooling 
timescale. In this paper we show that AGN winds can shock 
heat gas in the ISM to above the escape speed - this con- 
tributes to driving a galactic wind. To ensure that the cool- 
ing of the shock heated gas is correct, we modified the ISM 
model of SH03 to better match the expected cooling rate for 
gas with roughly a solar metallicity. To compute the local 



cooling rate we use a fit to Sutherland & Dopita ( 1993 1 for 



gas with solar metallicity. The fit is essentially eq. (12) of 
[Sharma et al.| ( |2010[ ) but with the cooling rate increased by 
a factor of two when the temperature is between ~ 3 x 10* K 
and 3 x 10^ K (to account for solar metallicity). From this 
cooling rate we compute a cooling time, tcooi', differences 
between the thermal energy and the subgrid "Ue// decay on 
treiax = tcooi- FoT Hiost of the relevant range of temperature 
and density in our simulations, this cooling time is shorter 
than the relaxation time of SH03. 

In addition to the above modification to the cooling 
rate, we also consider the role of inverse Compton cool- 
ing/heating. To do so, we modify the cooling/relaxation 
timescale of the gas to be 



t; 



— ^rnnl + ^C ' 



where the (non-relativistic) Compton time is given by 



tc = lOVr 



f R 



V300pc 



L 



10*6 erg s- 



(7) 



(8) 



with R the distance to the BH and L the associated AGN lu- 
minosity. We do not consider any radiative transfer effects in 
this paper, which in reality can modify both the luminosity 
L seen at a given radius and the local Compton temperature 
at a given radius. Instead, we take the Compton temperature 
to be that appropriate for the mean spectrum of luminous 
AGN, including the effects of obscuration: To ^ 2 x 10^ K 



(Sazonov et al. 20041 



In the limit that the Compton timescale is short com- 
pared to the two-body cooling time, we no longer relax the 
thermal energy of the gas to the effective equation of state 
value for the energy, Ue//, but rather to uc = ?)/2kTc- To 



transition between these two limits, we in general let the 
thermal energy of gas relax to: 



^relax 



tcoolUC + tgUeff 
tc + tcooi 



(9) 



For dense gas, atomic cooling dominates and the gas 
rather quickly approaches the sound speed associated with 
the effective equation of state. For gas densities character- 
istic of the ISM in the central kpc of our model galaxies 
(~ 10 — 10'^ cm~'^) the sound speeds are ~ 40kms~^ and 
thus the gas is primarily rotationally supported rather than 
pressure supported. This justifies our use of the viscous ac- 
cretion rate in equation [T] 



3 RESULTS 

Table IT] lists the parameters of the galaxy merger simula- 
tions presented in this work. We consider a single galaxy 
mass and merger orbit (see 



2.1 1, and explore the effects of 



including AGN wind feedback, both with and without radi- 
ation pressure. For the AGN winds, we vary both the total 
momentum flux in the wind (r^,) and the wind speed v-u,, 
with fiducial values of Tu, = 5 and v-w = 10, 000kms~^. 

In varying the AGN wind parameters, we are not guar- 
anteed that the resulting BH mass will be consistent with 
the Mbh — o" relation. In such cases, the effects of the AGN 
wind might not be realistic because of the unphysical BH 
mass. To quantify this, we present results with and without 
the radiation pressure feedback model explored in DQM. 
Our simple model of radiation pressure feedback produces 
model galaxies roughly on the observed Mbh — cr correlation 
for r ~ 20, the value used here (DQM). For the present pur- 
poses these calculations are useful primarily because they 
allow us to study the effects of AGN wind feedback for sys- 
tems in which the BH is guaranteed to be approximately on 
the Mbh — o" relation. We also separately carry out simula- 
tions with wind feedback alone. 



3.1 Effects of AGN Winds 

3.1.1 Black hole accretion and star formation rates 

Figure IT] shows the net accretion rate into the black hole, 
Min, and the star formation rate, M*, as a function of 
time (summed over both galaxies) for three different sim- 
ulations: the fiducial simulation (black curve) with radia- 
tion pressure (t = 20) and AGN wind feedback (t^ = 5 
and Vw = 10, OOOkms"^), a run identical to the fiducial run 
but without the wind feedback (gray), and a run with both 
wind and radiation pressure feedback that includes Compton 
heating/cooling in the thermodynamics of the gas (blue). We 
first describe the effects of the AGN wind and later return 
to the role of Compton heating/cooling. 

The inclusion of the wind feedback has little effect on 
either the accretion or star formation rate before the final co- 
alescence of the two galaxies at f ~ 1.7 Gyr. This is a priori 
surprising because equation (Im implies that for a given set 
of conditions at large radii (that determines Mvisc) the BH 
accretion rate at small radii is a factor of 1 + T^ric/v^ — 16 
smaller for the simulation with AGN winds than for the run 
with just radiation pressure. After a small number of time 
steps, however, the feedback due to radiation pressure is so 
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Table 1. Simulation Parameters 



Run 

Name 


Tw 


[kms-1] 


ISM 
Model 


^BHJ 
[IO^Mq] 


[10»Mq] 


[IO^Mq] 


/disk 


/out 


Mout 

[IO^Mq] 


Notes 


wol-10 


1 


10000 


a 


56.5 


8.67 


1.79 


0.014 


0.343 


5.47 


wo = wind only 


WO3-10 


3 


10000 


a 


17.9 


7.35 


1.05 


0.002 


0.472 


7.52 




woS-lO 


5 


10000 


a 


12.2 


7.36 


1.02 


0.002 


0.477 


7.60 




wo5-7 


5 


7000 


a 


5.36 


10.7 


2.91 


0.041 


0.175 


2.78 




wo5-3 


5 


3000 


a 


6.65 


12.1 


5.23 


0.059 


0.056 


0.89 




wolO-10 


10 


10000 


a 


6.20 


7.35 


0.92 


0.002 


0.466 


7.42 




wolO-7 


10 


7000 


a 


2.99 


10.8 


2.55 


0.014 


0.191 


3.04 




wo 10-3 


10 


3000 


a 


3.96 


12.2 


4.70 


0.069 


0.060 


0.96 




wp5-10 


5 


10000 


a 


4.38 


10.9 


3.11 


0.006 


0.196 


3.13 


fiducial sim; wp = wind + prad 


wp5-7 


5 


7000 


a 


5.04 


10.9 


3.77 


0.016 


0.175 


2.79 




wp5-3 


5 


3000 


a 


2.63 


12.8 


4.02 


0.159 


0.102 


0.36 




wpl-10 


1 


10000 


a 


11.8 


12.4 


3.86 


0.006 


0.100 


1.61 




wp3-10 


3 


10000 


a 


6.15 


11.9 


4.30 


0.006 


0.149 


2.37 




wplO-10 


10 


10000 


a 


3.36 


10.3 


3.21 


0.002 


0.268 


4.26 




wpfgSO 


5 


10000 


a 


10.2 


39.5 


4.54 


0.003 


0.148 


7.06 


/9,« = 0.3 


wpsli 


5 


10000 


c 


6.14 


10.8 


3.41 


0.002 


0.234 


3.73 




wpcomp 


5 


10000 


b 


6.30 


10.8 


3.21 


0.004 


0.276 


4.40 




wpcompl 


5 


10000 


b* 


5.01 


11.4 


3.96 


0.007 


0.181 


2.88 




wpinst 


5 


10000 


d 


5.90 


10.4 


2.67 


0.018 


0.216 


3.44 




wplr 


5 


10000 


a 


5.28 


10.2 


4.45 


0.004 


0.307 


4.89 


Lower resolution 


wplir 


5 


10000 


a 


4.77 


11.5 


3.94 


0.012 


0.166 


2.65 


Higher resolution 


po 








a 


9.33 


13.6 


4.92 


0.051 


0.007 


0.11 


No wind feedback 


pofg30 








a 


25.9 


44.5 


6.45 


0.028 


0.007 


0.32 


fg,i = 0.3, No wind feedback 



Model parameters and key quantities in our galaxy merger simulations. The top section lists runs with only the AGN wind feedback 

model (parametrized by the wind speed Dm and momentum flux t^). The bottom section lists runs with both the AGN wind model 

and t he radiation pressure model with r = 20. ISM model refers to one of four following options: (a) our fiducial cooling model (see 

j |2.5| , (b) same as (a) but including Compton heating with Tq = 2 x 10^ K, (b*) same as (a) but including Compton heating with 

Tc = W^ K, (c) same as in SH03, and (d) the gas returns to the effective EOS of SH03 instantaneously. Other columns list the key 

properties of the black holes, stars, and gas in the simulations: Mbh.j is the final black hole mass, M^f is the total mass of new stars 

formed by the end of the simulation, Mt^p is the total mass of new stars formed after the peak of star formation, /disk is the fraction of 

the gas mass within 3 kpc of the BH at the end of the simulation, /out is the fraction of the gas mass at large distances (> 10 kpc) 

above the orbital plane at the end of the simulation, and Afout is the corresponding gas mass, /disk is a proxy for how much gas 

remains in the galaxy at late times while /out is a proxy for the fraction of the gas that is in the outflow at late times; both /disk a-^d 

/out are normalized to the initial total gas mass of 1.59 X IO^^Mq. 



effective in all of the simulations in Figure [T] that the phys- 
ical conditions at small radii quickly adjust so that there 
is a balance between radiation pressure and gravity. This 
sets the BH accretion rate Mi„ to be ~ fgO^KTrjcG) (De- 
Buhr et al.||2010 1 independent of the presence of the AGN 
wind removing mass from the nuclear region (where fg and 
a are the gas fraction and velocity dispersion at small radii, 
respectively) . 

Although the AGN wind feedback has little effect on the 
early time star formation and BH accretion. Figure [T] shows 
that after final coalescence both the BH accretion rate and 
the star formation rate decrease significantly more rapidly 
with the inclusion of the AGN wind. The effects of the BHs 
feedback are maximized at and after final coalescence be- 
cause this is when the BH reaches its final mass (to within 
a factor of a few) and when the BH accretion rate is largest. 
Thus, the dynamics at and after final coalescence is the most 
sensitive to the details of the feedback physics. 

To assess why the late-time accretion and star forma- 
tion are suppressed by the inclusion of the AGN wind. Fig- 
ure [2] shows the gas mass fraction (normalized to the gas 
mass at the start of the simulation) within 3 kpc of the 
BH as a function of time, for the same three simulations 



as in Figure Ip Tidal interactions drive gas into the center 
of the two galaxies (or the merged remnant) after the first 
close passage (f ~ 0.5 Gyr) and at final coalescence (t ~ 1.7 
Gyr), maintaining a large nuclear gas fraction in spite of the 
intense star formation. The resulting increase in the BH ac- 
cretion rate (seen also in Fig. [I]) significantly increases the 
strength of the AGN feedback. In the simulation including 
AGN winds, this feedback efficiently removes material from 
the central regions of the galaxy, reducing the amount of gas 
inside the central 3 kpc of the merged system by a factor of 
~ 10 relative to the simulation without explicit AGN winds. 
Figure [3] shows images of the projected gas density for 
the fiducial simulation (/e/i) and the corresponding simula- 
tion without the AGN wind [right) &t t = 1.71 Gyr, shortly 
after final coalescence. The images are 280 kpc on a side 
and show a roughly edge-on view of the orbital plane of the 
two galaxies. While there is material at large distances in 
the simulation without the AGN wind, these features are 
the tidal tails generated during the merger and most of this 



^ The precise choice of 3 kpc is somewhat arbitrary but the re- 
sults in Figure [2] are relatively robust to changes in this choice. 
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Figure 1. The star formation rate, M*, and accretion rate into 
the black hole, Mi„, as a function of time for the fiducial simula- 
tion with both AGN wind and radiation pressure feedback (run 
'wp5-10' in Table 1; black), the run with only the radiation pres- 
sure feedback (run 'po'; gray), and the run that included Compton 
heating (run 'wpcomp'; blue). 
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Figure 2. The fraction of the total gas mass inside a sphere of 
radius 3 kpc (relative to the total initial gas mass) centered on 
the black hole for the same three runs shown in Fig. 1. When 
the AGN wind model is included, the central material is rapidly 
ejected after final coalescence (t ^ 1.7 Gyr). 



material is near the orbital plane. Thus our simple imple- 
mentation of radiation pressure feedback is not efficient at 
unbinding gas from the galaxy, though it is very effective at 
regulating the growth of the BH itself (DQM). By contrast, 
in the simulation with the AGN wind. Figure Is] (left panel) 
shows that there is significantly more material blown out of 
the galaxy, especially in the directions perpendicular to the 
orbital plane. Quantitatively, at the end of the simulation 
the mass of gas at large distances {\z\ > 10 kpc) from the 
orbital plane is ~ 3 x lO^M© in the simulation including 
the AGN wind, about 20 times larger than in the simula- 
tion with only radiation pressure feedback (see Table [T] and 




Figure 3. The projected gas density for the fiducial simulation 
(left) and the run without AGN wind feedback (right) at a time 
just after the final merger (t = 1.71 Gyr). Brighter color corre- 
sponds to higher density. The images are edge-on to the plane 
of the orbit and the box size is 280 kpc. While both simulations 
have material at large distance in tidal tails, there is significantly 
more material out of the orbital plane with AGN wind feedback. 
These images were generated using SPLASH ||Price|2007[l. 
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Figure 4. Effects of AGN winds on the distribution of radial 
velocities of the gas particles (relative to the black hole) at the 
end of the simulation. Few gas particles have velocities greater 
than ~ 200kms~^ in the run with only radiation pressure (run 
'po'; grey). By contrast, in the fiducial simulation with the AGN 
wind, a significant fraction of the particles have velocities above 
a few 100 kms~^. About 25% of these outflowing particles were 
explicitly added to the wind and kicked with an initial wind ve- 
locity of 10,000 kms~^ (blue), while the rest were not explicitly 
added to the wind (black) initially but were accelerated due to 
hydrodynamic interactions with the 'kicked' high-speed particles. 



Fig. [m below). A corollary of this efficient removal of gas by 
the AGN wind is that the total stellar mass formed during 
the simulation is ~ 20% smaller in the case with the AGN 
wind; most of this suppression in star formation happens at 
late times, after the final coalescence of the two galaxies (see 

Fig.g. 



3.1.2 Impact on the ISM 

The evacuation of the central part of the galaxy by the AGN 
wind would not be surprising if the majority of the material 
ejected was explicitly added to the wind by our deposition of 
momentum (i.e., if the unbound mass was primarily mate- 
rial that was explicitly 'kicked'). We find that this is not the 
case. At the end of the fiducial simulation (t = 2.85 Gyr), of 
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Figure 5. Effects of AGN winds on the mass weighted distribu- 
tion of gas temperatures. The same two simulations as in Fig. 4 
are shown at two times: near the peak of accretion at 1.71 Gyr 
(top), and at the end of the simulation at 2.85 Gyr {bottom). The 
temperature distributions are normalized to have a unit area. The 
excess of gas at high temperatures near the peak of BH accretion 
{top) in the case with wind feedback is due to the shock heating 
of the ISM gas. 



the 33,028 gas particles at large transverse position (l^l > 10 
kpc), only 8,904 of the particles have been explicitly kicked. 
The other 72% (82% by mass) have been ejected because 
of hydrodynamic interactions with wind material. In addi- 
tion, at any time only a small fraction ~ 5% of the gas 
inside Race (the accretion/feedback region) has been explic- 
itly added to the wind. These results demonstrate that the 
majority of the galactic outflow arises self-consistently due 
to the interaction between the AGN wind that we initialize 
and the surrounding ISM. 

Figure |4] illustrates this point in a different way: here 
we show the radial velocities (relative to the BH) of the 
gas particles at the end of the simulation, without (gray) 
and with (black and blue) wind feedback. Those parti- 
cles that have not been explicitly kicked and added to the 
wind are shown in black while those that have been explic- 
itly kicked are shown in blue. In the case with the AGN 
wind, there is a significant fraction of the gas with velocities 
~ 300 — 1000 km s~^ , larger than any of the gas in the simu- 
lation without the wind feedback (gray). Moreover, the out- 
flow velocities are similar for particles that have, and have 
not, been explicitly kicked. The particles that are explicitly 
kicked receive initial impulses of «„, = 10, 000 km s^^ and 
yet have final speeds of a few thousand kms~^. This result 
indicates that the wind we initialize sweeps up a significant 
amount of material in the ISM of the surrounding galaxy, 



driving additional gas out of the galaxy at high velocities 
~ 300 -1000 km s~^ 

Figure [5] shows the mass weighted temperature distri- 
bution (normalized to the total gas mass) near the peak of 
accretion (f = 1.71 Gyr; top) and at the end of the simu- 
lation (t = 2.85 Gyr; bottom) for the run with AGN wind 
feedback (black) and without (gray). The presence of the 
AGN wind leads to a significant amount of mass heated to 
above ~ 10^ K when the BH is near its peak accretion rate 
{top panel in Fig.lsl). This is due to gas that has been shock 
heated by the AGN wind. Much of this shock heated gas is 
ultimately able to escape the host galaxy - this gas expands 
outwards and cools adiabatically, leading to the large excess 
of ~ 10^~® K gas at the end of the simulation with the AGN 
wind {bottom). 

In our calculations that include a simple treatment of 
inverse Compton scattering off of the AGN's radiation field, 
we take the Compton temperature to be Tc — 2x 10^ K, the 
average observed value ( Sazonov et al.|2004 1 . This is similar 
to the peak in the temperature distribution we find in sim- 
ulations without Compton scattering {top panel of Fig. [5|. 
As a result Compton scattering mildly suppresses the high 
temperature tail of the gas shock heated by the AGN wind. 
This in turn suppresses the amount of mass blown out of 
the inner regions of the galaxy at late times (see Fig. [2|. 
Nonetheless, there is still a factor of ~ 3 — 10 times less gas 
in the inner 3 kpc in the simulation with the AGN wind 
and Compton cooling relative to the simulation without the 
AGN wind. In addition, the star formation rate and BH ac- 
cretion rate are significantly lower at late times even when 
Compton cooling is taken into account (Fig. fTJ . 



3.2 Dependence on Feedback Model Parameters 

In the previous section we demonstrated that an AGN wind 
generated at small radii can produce a galaxy-wide outflow 
that unbinds several 10^ Mq of gas and suppresses star for- 
mation in its host galaxy. In this section we explore in more 
detail how this phenomenon depends on the properties of 
the AGN wind. 

To quantify the effects of the AGN wind, we calculate 
the following quantities in all of our simulations and list the 
results in Table [l] (1) the final BH mass Mbhj, (2) the 
total mass of new stars formed by the end of the simulation, 
i\f,./, (3) the total mass of new stars formed after the peak 
of star formationr] M,,p, (4) the fraction of the gas mass 
within 3 kpc of the BH at the end of the simulation, /disk ~ 
Mg^rOkpc/Mg^tot, and (5) the fraction of the gas mass at 
large distances (> 10 kpc) above the orbital plane at the end 

The fractions 



3, Iz^lOkpc/JMs, tot 



/Mg 



of the simulation, /out = M, 
/disk and /out are both normalized to the initial total gas 
mass of 1.59 x W^'^Mq (from both galaxies); /disk is a proxy 
for how much gas remains in the galaxy while /out is a proxy 
for the fraction of the gas that is in the outflow at late times. 
Table [l] also quotes Mout, the absolute mass in the outflow 
at the end of the simulation. 

Figure [6] compares these key quantities for sinmlations 



^ We distinguish between M»_p and M^j because Figurefllshows 
that the AGN wind's largest effect on the star formation history 
is at late times, after the peak of BH accretion and star formation. 
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utilizing only the AGN wind feedback model (triangles) and 
for simulations with both radiation pressure and AGN wind 
(crosses). For comparison, we also show the results of the 
simulation with only radiation pressure feedback, taking r = 
20 (dotted lines). In the left set of panels, we fix t^ = 5 and 
vary the wind speed from 3, 000 — 10, OOOkms"^. In the right 
set of panels, we fix v^ ~ 10, 000 km s"'^ and vary t^ from 
1-10. 

Figure [6] {left panel) shows that the final BH mass is 
relatively independent of the speed of the AGN wind for 
a given value of r„,. By contrast, for simulations with the 
AGN wind model alone, the right panel shows that the final 
BH mass decreases by a factor of ~ 10 as t^ increases from 
1 — lOljThis is consistent with the analytic results of King 



( 2003[ ) and Murray et al. (20051, and the numerical results 
of DQM, in which the final BH mass decreases linearly with 
the total momentum supplied by the AGN, be it in the form 
of a wind or radiation. In particular, when the accretion 
rate is Eddington limited, the net force associated with the 
feedback is oc r^L^dd oc t^Mbh (eq. H|. Thus the BH mass 
at which feedback is able to overcome the gravity of the gas 
in the galactic nuclei scales as Mbh oc t~^ . 

Although the BH mass is only a weak function of «„ 
at fixed Tm, the properties of the gas at the end of the 
simulation depend strongly on v^n (bottom, left panels in 
Fig.|6|. Specifically, a larger wind speed v^ at fixed r^ leads 
to more efficient blow out of gas from the galaxy at late 
times (smaller /disk and larger /out). For a given value of 
T„, the momentum fiux in the outflow (p„ = Tn,L/c; eq. |4]) 
is independent of «„ . By contrast, the energy fiux in the out- 
fiow {Ew — pwVw) increases for a larger wind speed. Thus a 
larger wind speed leads to more efficient shock heating of the 
gas in the galactic nucleus. A larger value of r^, has a simi- 
lar effect, increasing both the momentum and energy fiuxes 
in the wind, and producing a more powerful galaxy-scale 
outflow (bottom, right panels in Fig. Im. 

Overall, then, the galaxy-wide outflow is a consequence 
of both the wind launched from small radii sweeping up and 
driving out ambient gas and the wind generating a shock 
that heats and unbinds ambient gas. This is why the gas 
content in the outflow - as measured by /out in Figure [6] 
- increases with both increasing u,„ and increasing t^. In 
general the simulations with AGN wind feedback alone have 
somewhat higher outflow rates than the simulations with 
AGN winds and radiation pressure. The reason is that the 
feedback is then provided entirely by the wind, enabling it 
to be more effective and unbind a larger amount of gas. 

The mass of new stars formed in our simulations de- 
creases for larger Uu, and/or t^, (top panels in Fig. p|. This 
is a simple consequence of the galaxy-scale outflow being 
more powerful for larger v^ and/or t^. The net effect of 
the outflow on the stellar mass formed in the simulation is, 
however, relatively modest (typically ~ 10-20%). The AGN 
wind's most prominent effect on the star formation history 
is that it can suppress the late-time star formation, as high- 
lighted in Figure [l] It is also important to note, however, 
that the outflow does not always have this effect. For ex- 



^ The same effect is much less evident in the simulations with 
radiation pressure and the AGN wind (crosses in the right panels) 
because the former largely sets the BH mass. 



ample, the top left panel of Figure [6] shows that simulations 
with outflows having Tm = 5 and v^ = 3000 km s~^ form a 
very similar amount of new stars (both in total and after 
the peak of star formation) to simulations with the radia- 
tion pressure feedback model alone. This is consistent with 
the fact that for these parameters, the outflow does not ap- 
preciably change the gas content in the galactic disk (/disk 
in Fig. Im), although it does unbind a modest amount of ad- 
ditional gas ~ 3 X IO^Mq (see /out in Fig. l6| and Mout in 
Table [l]). 

3.3 Numerical Tests 



We argued in j ]3.1| and |3.2| that the galactic outflow seen 
in our simulations is a consequence of both the AGN wind 
sweeping up ambient gas (momentum conservation) and the 
AGN wind shock heating ambient gas and unbinding it (en- 
ergy conservation). The former mechanism does not depend 
sensitively on the thermodynamics of the gas, while the lat- 
ter does. As described in §2.5| we modified the ISM equation 
of state of SH03 in order to include two-body cooling at solar 
metallicity; we also included a simple treatment of inverse 
Compton scattering off of the AGN's radiation field in some 
of our simulations (see Fig. [I]). 

To further assess the sensitivity of our results to the 
thermodynamics of the gas, we carried out a number of addi- 
tional tests. In particular, two extremes are to use the SH03 
model as is (in which gas that deviates from the equilibrium 
equation of state relaxes back on a timescale different from, 
and typically longer than, the true cooling time), or to force 
gas to return to the effective equation of state of SH03 in- 
stantaneously (thus removing the possibility of strong shock 
heating unbinding gas). The results of these calculations 
with our standard AGN wind parameters are given in Ta- 
ble IT] labeled wpsh and wpinst, respectively. In both cases, 
the AGN wind drives a galaxy-scale outflow with an ejected 
mass ~ 3 X IO^Mq, very similar to our fiducial simulation. 
However, these 'extreme' ISM models have effects that are 
not surprising: SH03's equation of state (in which the gas ef- 
fectively cools more slowly) is more effective at suppressing 
star formation and clearing gas out of the nucleus (wpsh), 
while the simulation with instantaneous cooling (wpinst) re- 
tains a larger nuclear gas disk at the end. Overall, however, 
these differences are modest - they are smaller, e.g., than the 
uncertainties introduced by the relatively poor current con- 
straints on the total nromentum/energy flux in AGN winds 
(tu,). This test demonstrates that our treatnrent of the ISM 
thermodynamics is probably sufficiently accurate to assess 
the impact of AGN winds at the order of magnitude level. 
However, future calculations that incorporate a more real- 
istic multi-phase, turbulent ISM (e.g., |Hopkins et al.||2011 l 
will be very valuable. 



Recently, Durier & Dalla Vecchia (20111 showed that 



after injection of energy, either thermally or kinetically, in- 
accuracies can arise in the subsequent evolution of the SPH 
particles when each is assigned a separate timestep, as oc- 
curs in GADGET. To ascertain if our feedback is subject to 
the innacuracies they report, we reran the fiducial simula- 
tion with an implementation of their timestep limiter. All 
of our key qualitative conclusions are unchanged, though 
there is a modest (~ 30 %) change in the post merger star 
formation and disk fraction (/disk). 
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Figure 6. Key quantities characterizing the effects of the AGN wind model as a function of wind speed v^ (left column) and momentum 
flux Tm (right column); see Table 1 for more details. These results quantify describe how the properties of the galaxy merger remnant 
at the end of the simulation depend on the parameters characterizing the AGN wind at small radii. In each panel, the triangles label 
the runs with only the AGN wind feedback, the crosses label the runs with both AGN wind and radiation pressure feedback, and the 
dashed line is for the simulation with only radiation pressure feedback. The four rows (from top to down) show: the mass Mt,p of new 
stars formed after the peak of star formation during the simulation; the final black hole mass Mguj; the mass fraction /disk of gas at 
the end of the simulation within 3 kpc of the black hole, and the gas mass fraction /out at the end of the simulation that is at large 
distances from the orbital plane: \z\ > 10 kpc . Both /disk and /out are normalized to the total initial gas mass in the system. 



Figure It] shows the radial velocities of the SPH particles 
(relative to the BH) at the end of the simulation for three 
runs that differ only in the mass resolution with 1.6, 4.8, and 
16.0 X 10^ particles, respectively. For the two higher resolu- 
tion simulations, the distributions are quite similar for the 
bulkofthedisk(«~ 0) and outflow (u ~ 500-1000 kms'^), 
though the high velocity tails differ somewhat. Table [2] 
compares the integrated quantities shown in Figure |6] The 
biggest difference is that the highest resolution simulation 
retains more mass in a nuclear disk at late times. However, 



the mass in the outflow and the fraction of the mass in the 
outflow at late times are reasonably converged. The total 
mass of new stars formed and the final BH mass also do 
not change appreciably with resolution. We conclude that 
the effects of this model are reasonably well resolved at the 
fiducial resolution employed. 
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Figure 7. The distribution of radial velocities relative to the 
black hole for the gas particles in the fiducial simulation (black), 
the simulation at lower resolution (red) and the simulation at 
higher resolution (blue). 



Table 2. Resolution Test 



Property 


High 


Fid 


Low 


Af,,t [IO^Mq] 


11.5 


10.9 


10.2 


Mbhj [IO^Mq] 


4.77 


4.38 


5.28 


/disk 


0.012 


0.006 


0.004 


/out 


0.166 


0.196 


0.307 


A/out [IO^Mq] 


2.65 


3.13 


4.89 



Simulation properties are the same as those of the fiducial sim- 
ulation in Table 1. The three simulations differ only in the mass 
resolution, with a total of 1.6 X 10^ (low), 4.8 X 10^ (fid), and 
1.6 X 10^ (high) particles. 

4 DISCUSSION AND CONCLUSIONS 

We have carried out three-dimensional SPH simulations of 
the interaction between a high-speed outflow produced by an 
AGN and interstellar gas in the AGN's host galaxy. We show 
that this interaction can drive a large-scale galactic outflow 
that in some cases unbinds a signiflcant fraction of the ISM 
of the host galaxy (see, e.g.. King 2005 for closely related ar- 
guments). The AGN-driven galactic winds found here pro- 
vide a possible explanation for the high velocity outflows 



observed in some post-starburst galaxies (Tremonti et al 



20071 and for the massive atomic and molecular outflows 
^ seen in local ultra-luminous infrared 



with V r^l, 000 km s 

galaxies (ULIRGs) ([Feruglio et al.]|2010[ iChung et al.|20rT 



Rupke fc Veilleux 2011 Sturm et al. 2011 1 . We return to this 



below. Our specific calculations assume that AGN activity 
is triggered by major mergers of galaxies (as seen, for exam- 



ple, by Koss et al. (2010l), but we suspect that the results 
presented here are much more general and apply relatively 
independent of the physical mechanism(s) that drive gas into 
the central ~ 100 pc of galaxies. The physically-motivated 
model utilized here for AGN wind feedback is easy to imple- 



ment (i2.4| and could readily be applied in other contexts. 



including high resolution simulations of galactic nuclei and 
cosmological simulations of BH growth and evolution. 

Our results provide a quantitative mapping between 
the properties of an AGN wind in the vicinity of the BH 
and the resulting large-scale galaxy-wide outflow (Table 
[I] and Fig. Im). We have parametrized the AGN wind at 



small radii in terms of its speed Um and momentum flux 
Pw ~ T^L/c - the corresponding energy flux in the wind is 
Ew = 0.5 p™ v^ = 0.5r,i,(u„/c)L. 

We find that an AGN wind with u„ ~ 10, 000 km s"^ - 
typical of, e.g., broad-absorption line quasars - produces a 
galactic outflow with velocities ~ 1, OOOkms"^. This reduc- 
tion in velocity is a consequence of the AGN wind sweeping 
up and driving out ~ 5 — 10 times as much mass as was in 
the wind initially (e.g.. Fig. p|. Quantitatively, we find that 
the AGN-driven galactic outflow unbinds ~ 1 — 5 x 10^A/q of 
gas, ~ 10 — 40% of the initial gas in the two merging galax- 
ies (which are Milky Way- like spirals in our calculations). 
This is true for a wide range of outflow models, covering 
r„ ~ 1 - 10, «„ ~ 7,000 - 10,000 kms"^ and with and 
without the radiation pressure feedback model from |DeBuhr| 
et al.| ( [20Tol|2011[ ). Not surprisingly, AGN winds with larger 
momentum and energy fluxes (larger t-u, and v^) are more 
effective at driving galaxy-scale outflows (Fig.™. 

The mass outflow rate from the vicinity of the galaxy 
is somewhat noisy and thus difficult to directly measure in 
the simulations. The average outflow rate relative to the 
star formation rate can, however, be readily estimated by 
comparing the total unbound mass (Mout in Table [I| with 
the mass in new stars formed in the simulation (M*j and 
Mt,p in Table [l] where M,j is the total mass of new stars 
formed and Mt,p is the mass of new stars formed after the 
peak of star formation). This comparison implies that the 
average mass outflow rate is typically comparable to the 
star formation rate; after the peak of star formation in the 
merger, which is when the BH reaches its final mass, the 
average outfiow rate often exceeds the corresponding star 
formation rate by factors of a few. 

Physically, we find that the properties of the galaxy- 
wide outflow are determined by both the momentum and 
energy fluxes in the AGN wind at small radii. A larger mo- 
mentum flux implies that more mass can be swept out of the 
galaxy while a larger energy flux in the wind (e.g., larger v^ 
at fixed pm) leads to more shock heating of the ambient ISM. 
This shocked gas partially cools but retains enough thermal 
energy to contribute to driving the galaxy-wide outflow. For 
quantitatively accurate results we find that it is necessary 
to modify the Springel & Hernquist (20031) equation of state 



to properly include metal-line cooling and Compton heat- 
ing/cooling as processes by which gas relaxes back to the 
effective equation of state that is used to calculate the ISM 
pressure (see j ]2.5[ |. 

To produce a galaxy at the end of the simulation that 
is on the Mbh — o" relation using only feedback by our 
AGN wind model, we require Tm > 5 — 10 (see Table [l] 
and Fig. [6] the velocity dispersion of the remnant galaxy 
is ~ 170kms~^ so that Mbh ~ 7 x IO^'M© for galax- 
ies on the mean Mbh — o" relation). This required mo- 
mentum fiux is comparable to, although somewhat smaller 
than, the required momentum flux of p ~ 20 L/c that DQM 
found using a simple model of radiation pressure feedback 
in which all particles in the vicinity of the BH feel a force 
oc r L/c. The energy flux corresponding to r„ ~ 5 — 10 and 
v-w ~ 10, 000 km s~^ is £^ ~ 0.1 L, similar to the energy in- 
jection rate of ~ 0.05 L that Di Matteo et al. (20051 found 



was required to reproduce the Mbh — o" relation (although 
they deposited thermal energy into ambient gas while we 
supply kinetic energy to a wind). 
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One important distinction that we find is between AGN 
driving galactic outflows and AGN outflows regulating the 
growth of the BH. For example, in our models with r„ = 1 
and Vm — 10, 000 km s~^, the flnal BH mass is a factor of 
~ 10 above the mean Mbh — o" relation in calculations that 
only include AGN wind feedback. Also including the simple 
radiation pressure feedback model from DQM in the simu- 
lation resolves this discrepancy; the same AGN wind model 
then unbinds ~ 10^ M© of gas, ~ 10% of the initial gas 
mass in the system. However, in spite of this outflow, the 
AGN wind neither regulates the BH growth nor has a sig- 
niflcant effect on the average star formation in the system. 
This demonstrates that very powerful AGN outflows with 
Tu, ~ 3 — 10 are required for the BH to have a substantial 
effect on either its own accretion history or the star forma- 
tion in the surrounding galaxy. This conclusion is similar to 



that of Silk & Nusser (20101, who required a comparable 



momentum to eject gas from their model galaxies, and to 



that of Kaviraj et al. (2011 1, who found that simple models 



for the color evolution of early type galaxies required feed- 
back similar in strength {E of a few percent) and duration 
(~ 200 Myr) to what we find here. 

The galactic outfiows found in our simulations with 
Vm ~ 10, 000 km s~^ and Tm > 3 have properties broadly 
similar to those observed in some massive post-starburst 
galaxies ([Tremonti et al. 20071 and in local ULIRGs (Fer 



uglio et al.|2010| [Chung et al.|2011||Rupke fc Veilleux|2011 



Sturm et al.|20lT |. Moreover, it is striking that Sturm et al.| 
(20111 inferred p^ ~ lOL/c, i.e., r™ ~ 10, for the molecu- 



lar outflows in local ULIRGs (these values are uncertain at 
the factor of ~ 3 level; see, e.g., [Sturm et alT||2011| for more 
details). This is comparable to the outflow momentum flux 
we flnd is required to both suppress late time star formation 
during galaxy mergers and produce remnants approximately 
on the observed Mbh — o relation. This correspondence is 
encouraging, though it does not directly address the ques- 
tion of what mechanism powers such outflows in the flrst 
place. The relative role of star formation and AGN in pow- 
ering the ~ 1000 km s"'^ molecular outflows is also unclear. 



with Sturm et al. (20111 presenting evidence from OH ob- 



servations in favor of AGN, while Chung et al. (2011 1 argue 



that star-formation is the dominant power source based on 
their CO observations. In addition, in some samples of post- 
starburst galaxies, the observed outflows appear consistent 
with those driven by star formation ( Coil et al.|2011 1. 

From the theoretical perspective, the biggest uncer- 
tainty in applying our results is the uncertainty in the ab- 
solute momentum/energy flux in AGN outflows at small 
radii. This uncertainty remains even for the most well- 
characterized outflows, those of broad-absorption line (BAL) 
quasars. In theoretical models in which BAL outflows are 
produced by radiation pressure on lines in the vicinity of 
the broad-line region, the predicted momentum flux is typ- 
ically < i/c, i.e., r^, < 1 ([Murray et"al]|1995[). This fol- 



lows from the fact that (observationally) the lines do not 
cover the entire continuum and thus cannot absorb the en- 
tire continuum momentum flux. This constraint is less clear 
for low-ionization BALs (in which the lines cover more of the 
continuum), but probably remains reasonably accurate. One 
way out of this conclusion is if magnetic flelds contribute to 
accelerating BAL outflows (as in, e.g.,iProga2003) since this 
can increase both the mass flux and terminal velocity of the 



outflow. There are observational estimates that the momen- 
tum flux in some low-ionization BAL outflows may indeed 
be > Ljc ( |Moe et a l. 2009, B autista et al.|2010|[Dunn et al 



20101, though these are difficult observations to interpret 



and they have complicated selection effects. 

In addition to line-driven winds, it is possible that 
even more powerful outfiows are driven during evolution- 
ary phases distinct from optically bright quasars. The two 
most promising from our point of view are: (1) outflows 
driven by radiation pressure on dust during phases when 
the ISM around the AGN is optically thick to far-infrared 
radiation ( [Murray et"aL][2005| [PeBuhr et aL]|2010| . In this 
case the momentum flux in a small-scale AGN wind can 
in principle reach ~ tL/c, where r is the infrared opti- 
cal depth of the galactic nucleus. Such an outflow is not 
directly seen in DQM's simple implementation of radiation 
pressure feedback. However, this could easily be a limitation 
of the lack of a proper treatment of the infrared radiative 
transfer. Indeed, we expect that radiation pressure will in- 
evitably drive a high-speed outflow off of the nuclear disk 
('torus') at ~ 1 — 10 pc. (2) An additional source of power- 
ful high-speed outflows can be produced if the fueling rate 
of the central AGN significantly exceeds Eddington (e.g., 
King|[2003 1 . In this case, the accretion disk becomes radia- 



tively inefficient and much of the nominally inflowing mass 



is likely unbound (e.g., Blandford & Begelman 19991. The 
disparity in spatial scales between the horizon and where 
the AGN accretion rate is set (~ 1 — 10 pc; e.g., [Hopkins] 
'& Quataert~2010) implies that it is very likely that there is 
a phase of super-Eddington fueling. The critical question is 
whether this phase lasts long enough to dominate the inte- 
grated properties of a black hole's outflow. 

All of the calculations presented in this paper utilize 
the effective equation of state model of [Springel fc Herii^ 
[quist| ( ^2003^ (modified as in j ^2.5[ ) to pressurize the ISM and 
prevent runaway gravitational fragmentation. This is a very 
simplistic treatment of ISM physics and in the future it is 
important to study the interaction between a central AGN 
and the surrounding ISM using more realistic models that 
approximate the turbulent, multi-phase structure of the ISM 
(e.g., [Dobbs et al.|20l"T| [Hopkins et al.|2011| ). 
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